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ABSTRACT 

Iron, the Universe's most abundant refractory element, is highly depleted in both circumstellar and 
interstellar environments, meaning it exists in solid form. The nature of this solid is unknown. In this 
Letter, we provide evidence that metallic iron grains are present around oxygen-rich AGB stars, where 
it is observationally manifest as a featureless mid-infrared excess. This identification is made using 
Spitzer Space Telescope observations of evolved globular cluster stars, where iron dust production 
appears ubiquitous and in some cases can be modelled as the only observed dust product. In this 
context, FeO is examined as the likely carrier for the 20-/im feature observed in some of these stars. 
Metallic iron appears to be an important part of the dust condensation sequence at low metallicity, and 
subsequently plays an influential role in the interstellar medium. We explore the stellar metallicities 
and luminosities at which iron formation is observed, and how the presence of iron affects the outflow 
and its chemistry. The conditions under which iron can provide sufficient opacity to drive a wind 
remain unclear. 
Subject headings: stars: mass-loss — circumstellar matter — infrared: stars — stars: winds, outflows 

— globular clusters: individual (NGC 362, NGC 5139, NGC 5927) — stars: AGB 

and post- AGB 



1. INTRODUCTION 

Iron is the sixth most abundant element in the Uni- 
verse, and the most abundant refractory element. It 
is observed to be highly depleted i n both interstel- 
lar and circumstellar environments ("Savage &: Bohlir 



200£ 



119791 : iMaas et all 120051 : [Felgado Inglada ct all 
and must therefore predominantly exist in an unknown 
solid form. Iron may be incorporated in other dust 
grains, pr imarily silicates, thou gh these are u sually 
iron-poor (jGail fc Sedlmavdll999HKemper et"alll2002f ). 
Alternatively, iron may form a metalli c condensate 
(|Kemper et al.ll200irVerhoelst et al.ll2009[ ). 

Dust grains form around evolved stars by condensa- 
tion from the gas phase, either directly, or onto molecu- 
lar 'seeds'. Individual dust species can usually be iden- 
tified by distinct infrared emission bands. In oxygen- 
rich environments, amorphous or crystalline forms of sil- 
icates, spinels and corundum are observed to form, whilst 
carbon-rich environments give rise to amorphous car- 
bon, graphite and SiC. Metallic iron grains, however, 
produce featureless infrared emission which can be dif- 
ficult to differentiate from other sources, particularly 
amorphous carbon dust. Iron has hitherto only been in- 
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ferred as a like l y component of oxygeii-rich dusty winds 
(|Kemper et all i2002t TVerhoelst et all 120091 ) . but never 

positively identified. 

Previous observ ational studies (jMcDonald et al.l 120091 : 
iBover et al.ll2009[ ) of giant stars in globular clusters have 
found a featureless contribution in addition to the flux 
emanating from the star's photosphere, which has been 
attributed to an unidentified circumstellar dust species. 
We herein show that this dust species is metallic iron. 



2. EVIDENCE FOR DUST EMISSION 
2.1. Observations 

Our sample consists of 35 highly- evolved, (strongly- 
)pulsating stars in globular clusters (jSloan et al.ll201Cl() 
(Paper I) , observed using the Spitzer Space Telescope In- 
frared Spectrograph (IRS). These long-period variable 
stars lie near or above the red giant branch tips of 19 
Galactic globular clusters, which range in metallicity 
from [Fe/H] = -1.62 to -0.10. Observed with Spitzer 's 
low-resolution modules, the spectra cover A = 5.2 — 38 
^m with a resolving power of i? '^ 60 — 120. 

One carbon-rich star (Lynga 7 VI) is discounted, leav- 
ing a sample of 34 stars. Of these, silicate features at 
10 and 20 /im are seen in 24 objects, implying oxygen- 
rich chemistry. More pertinently, the remaining ten stars 
either have very weak silicate features or appear to be 
'naked'. 

JHKL-ha.nd light curves are also available (Paper I) 
for 25 of our 34 objects, which allows us to construct 
the photometric fluxes of these objects at the pulsation 
phase of their IRS spectrum. 

Figure[T]shows the L-[8] colour for these objects, where 
the 8-fim photometry was derived by convolving the IRS 
spectrum with the Spitzer IRAC 8-fiin filter. A truly 
naked photosphere should have L — [8] w 0, once in- 
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Fig. 1. — Extinction-corrected L-[8] colours in our sources. Plus 
signs show stars with pure iron dust, filled circles show stars also 
hosting silicate dust (symbol size is proportional to strength of 
emission). Crosses show stars exhibiting a strong 20- /xm feature. 
Colour does not vary significantly with E{B — V), showing that 
this is not merely a reddening effect. 

terstellar reddeni ng; has been taken into accounU (c.f. 
iBover et al.ll2008[ ). All 25 stars, however, have a positive 
colour. This implies a flux excess at 8 /Ltm which can- 
not be caused by commonly-assumed dust species (sili- 
cates, alumina) under normal conditions of temperature 
and grain size. This excess flux at long wavelengths is 
confirmed by comparison to MARCS photosphcric model 
spect ra (|Gustafsson et all 119751 120081: iMcDonald et alj 
l2009f ) as shown below. 

From our sample, four example stars have been se- 
lected for detailed modelling: NGC 5927 V3, a highly- 
evolved AGB star with a strong, narrow, lO-fun sili- 
cate feature; NGC 6352 V5, a less luminous giant star 
with unusual dust features; and NGC 362 V2 and NGC 
5139 V42, 'naked' RGB-tip or AGB stars showing optical 
TiO bands (indicating they are oxygen -rich) and infrared 
dust excesses, but no silicate features (iSmith et al.l 119991 : 
IMcDonald et al.ll2009l: iBover et al.|[2009ir 

2.2. Fitting the observations 

For a dust-free comparison, we use a stellar atmo sphere 
model created using the MARCS code (Gustafsso n et all 
[T975lf 20081 at 3500 K, [F e/H] = -1.0, [a/Fe] = +0.3 
(see [McDonald et aLl 120091 for details of these models). 
Being very similar objects near the giant branch tip, we 
expect the effective surface temperature below the dust- 
producing zone to be 3750 ± ~250 K. We use the near- 
infrared flux to fit the photospheric continuum. A low- 
temperature, metal-rich model therefore represents an 
intentional bias toward a model with a more-luminous in- 
frared spectrum compared to its near-infrared flux. This 
reduces the calculated excess to rule out disparity be- 
tween the modelled and real photospheric spectra. 

Figure [2] shows the observed spectra, alongside 
fits to the stellar wind contribution using dusty 
(jNenkova et al.l [19991) . Using the MARCS model atmo- 
sphere as an input, we fit a stellar wind model us- 
ing metallic iron grains (|Ordal et al.l I1988D . For each 
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model we assume a radiatively-driven wind; and a stan- 
dard Mathis-Rump l-Nordsieck grain size distribution 
(jMathis et al.lll977[ ). given by n{a) = a~«, q = 3.5 for 
a = 5 — 250 nm. 

We correct for interstellar reddening using the absorp- 
tion profiles from McCl ure (200 9). We further assume 
that Ay = 3.2E{B - V) and that Ay/Ak = 7.75. 
EjB- V) va lues for NGC 362, 5927 and 6352 are from 
[Harris! ([199^ . E(B - V) = 0.10 is assumed for to Gen, 
based onlHarrid (|l996[ ) and IMcDonald et al.l ^QO^. The 
correction applied is thus quite small and negligibly dif- 
ferent from that applied in Paper I. 

NGC 362 V2 and u Cen V42 are best fit with a dust 
shell truncated to an outer radius (router) of '^100 times 
the inner radius (rjnnor = the dust formation radius). 
(Normally, fits are extended to router/nnnor }t 1000.) 
This suppresses the flux at A > 30//m. A similar ef- 
fect can also be achieved by using a grain distribution 
with fewer large grains: both effects may be in play (we 
only model a truncated shell here). The strength, and in- 
deed existence of, this suppression depends crucially on 
the subtraction of the infrared background (caused by 
Galactic interstellar medium and other cluster objects). 
While every effort has been made to make sure this has 
been done correctly, we cannot be conclusively sure that 
such a suppression and thus truncation exist. 

In NGC 5927 V3 and NGC 6352 V5, we also observe 
other dust species, which emit at A > 8/im. Iden- 
tified species include silicates (magnesium-rich olivines 
and pyroxenes) at 10 and 20 ^m and amorphous alu- 
mina (AI2O3 ) at 1 1~15 /im (optica l cons tants from 
iDraine fc Leg [19841: [Begemann et al.l I1997D . Related 
crystalline minerals with a smaller degree of amorphisa- 
tion are likely responsible for the sharp 9.6- and 13.1-/im 
features seen in NGC 6352 V5 (Sloan et al. 2003). Note 
that the fit to NGC 5927 V3 produces a fairly good match 
to the 10-/xm feature, but at a wavelength which is too 
red. This can be cause d by a number of factors, includiii g 
grain size and porosity ([Voshchinnikov fc: Hennin"^l2008l) . 
The 20-^m feature has an unknown carrier: a promis- 
ing candidat e is magncsio-wtistitc (Mg^jFci-^O; see, e.g., 
iPosch et al.l 120021 : [tebzcltcr et al. 2006: Paper I). This 
provides a good fit to the feature in NGC 6352 V5. The 
comparatively-red peak wavelength seen in NGC 6352 V5 
further suggests that a; ~ 0: i.e. that it is near-pure FeO 
([Henning et al.l[l995t ). We return to this feature later. 

We could not fit a single-temperature dust model to 
these two stars, implying that the dust species are not 
in thermal contact with each other. Iron and alumina 
were fit at temperatures of ~1000 K. In NGC 6352 V5, 
sihcates and FeO were fit at ~450 K; in NGC 5927 V3, 
the strong 10-/im peak suggests mainly glassy silicates, 
which we fit at --400 K. 

Table[T]lists our fitted parameters, including dust tem- 
peratures at the outer and inner edges of the dust enve- 
lope (Toutcr, Tinner) and the size of the envelope used 

(routcr/nnncr, Hnnor ~ 2 R* ) . Wc do not iucludc maSS- 

loss rates, as explained at the end of this section. 

2.3. Alternative sources of excess 

Featureless infrared excess can be produced by a va- 
riety of circumstellar material. This can take the form 
of free-free emission, emission from shells of molecular 
gas (a 'molsphere'; e.g. ,Tsuii.,200Qil . extremely-large sili- 
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Fig. 2. — Modelled dust contributions (solid lines) to Spitzer IRS spectra (black dots) and JHKL photometry (black pluses) which has 
been corrected to the pulsation phases of the IRS spectra. The dashed line is a MARCS model denoting the photospheric contribution. The 
bottom two panels include the extra indicated contributions from AI2O3, silicates and FeO (dotted lines). 



TABLE 1 

Parameters of dusty models presented in Figure [2] 



Parameter 



NGC 362 V2 u) Cen V42 NGC 5927 V3 



NGC 6352 V5 



Assumed luminosity (Lq) 1826 1862 2000 

Assumed [Fe/H] -1.16 -1.62 -0.37 

Assumed £(5 - V) 0.05 0.10 0.45 

References 1,3 2 4 



2000 
-0.70 
0.21 

4 



Wind component 



Composition 

-'inner dust envelope edge V / 
-'outer dust envelope edge K^J 
Pouter/ dinner 

Optical depth at 0.55 fiia 



100% 
iron 

800 

194 

50 

0.28 



100% 
iron 

1000 
173 
120 
0.70 



100% 
iron 

1100 

1000 
0.55 



100% 
silicate 

400 

1000 
0.45 



96% Fe 
4% AI2O3 

1000 

100 
1.05 



70% Silicate 
30% FeO 

450 

100 
0.12 



References: (ll lBover et"an iWO^ : (2}JMcDonald et al.l 120091 '): (3) E{B-V) from lHarrij irgg^l : (4) 

m [Harris [1996) . 



luminosity estimated, metallicity from[ 
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cate grains, or amorphous carbon grains. Figure [3] shows 
modeUed fits for these possible emission sources, mod- 
elled in an identical way to the iron fits listed above. 

lonisation in a chromosphere or shock fronts propa- 
gating through the outer atmosphere can give rise to 
free-free emission, as observed by radio ob servations to- 
ward Galactic Miras (|Reid fc Mentenlll997[ ). The optical 
depth of free-free emission in a region of thickness r can 
be calculated as: 



T^ == / 3.28x10" 



lio^Ky 



-1.35 



^GHz 



1 / »c \ 

' Vcm-3/ 



2 5r 
pc 



(1) 



where the electron density, rip, is recoverable from the 
electron temperature. To, via the Saha equation. At low 
frequencies, where Tjy ;» 1, the emission scales as F^ oc 



for an isothermal, isobaric region; and F^ oc v 



-0.6 

for an isothermal wind expanding with constant velocity. 
For uj Cen V42, we find our spectrum would require Ti, = 
1 at 3-4 /im, and that F^, oc z/~^-^ for t^, ;» 1. This more- 
closely approximates an isobaric density profile, which 
would be difficult to achieve in an expanding wind. Very 
similar conditions would have to exist in all 34 of our 
stars, as the shape of their spectra are very similar. 

An isobaric chromosphere producing a spectrum with 
this optical depth cannot exist except in the following 
case: the wind must be mostly ionised (> 6000 K) and 
either the chromosphere extends to an unphysical ra- 
dius (> 70 R0) and/or the mass- loss rate becomes un- 
physical (> 10~^'^ M0 yr^^) for a thermally-expanding 
wind. Departing from an isobaric state makes condi- 
tions even more prohibitive. While chromospheres may 
be present in these stars, this infrared excess is not seen 
in chromospherically-ac t ive stars that are less e volved 
([McDonald fc van LoonI 120071: iBover et all |200£i). In 
V42, circumstellar ionisation app ears d ominated by pul- 
sation shocks. Reid & Mentenj ()1997[ ) show that, for 
Galactic Miras, the temperature and density increases 
provided by shocks in the outer atmosphere are insuffi- 
cient to produce an ionised layer that would be optically 
thick in the infrared. We thus rule out free-free emission 
as a dominant contributor to our spectra. 

The molsphere creates a dense 'forest' of molecular 
lines. Water lines are the most important, and cause 
emission throughout the spectrum, but mainly in the 
6-8-/im and > 20/im regions. This emission is seen in 
some of our example stars (Paper I). However, the water 
bands in our spectra are very weak compared to star s 
with naked molspheres (e.g. /3 Peg — iTsuii et al.lll997t ). 
Gas-phase model ei nission spectra (Figure El) , c reated 
from the line list of [Partridge fc Schwenki (|1997[ ) using 
SpectraFactorjQ (|Cami et al.ll2010D , fail to reproduce the 
fiatness of the spectrum for any reasonable combination 
of temperature and column density. Figure |3] shows the 
highest column density we consider feasible, and even 
this contains ~ 10~^ Mq of water or ~0.01% of the 
star's oxygen. Furthermore, most molsphere host stars 
still require an extra opacity sou rce (such as iron grains) 
providing fiux at 6-8/ini (Vcrhoc lst et al.i r2009'). Infrared 
molecular emission should therefore not make a signifi- 
cant contribution to our spectra. 

Of the commonly-considered dust species only metal- 

* http://spectrafactory.net 



lie iron and amC produce a featureless infrared contin- 
uum that matches the observed 'naked' spectra. Other 
dust species (e.g. graphite) produce either strong spec- 
tral features or inflections which do not fit observed spec- 
tra. Spectral features can be damped by unusual condi- 
tions, notably exceptionally large grains (jII6fnerl i2008): 
these cannot simultaneously reproduce both the excess 
at A < 8^m and A > 20/xm at the observed levels (Fig- 
ure [3]) and are not expected in metal-poor stars where 
grain growth may be hampered by a lack of constituent 
elements. 

Obscrvationally, amC can also provide an acceptable 
fit to our example stars over most of the spectrum. The 
remaining discrepancies may be fixable by fine-tuning the 
wind parameters (grain size, density distributions, tem- 
peratures, background subtraction, etc.). Detection of 
amC as dominant dust species in an oxygen-rich envi- 
ronment would be unexpected. In giant stars, carbon 
and oxygen are bound into CO near the stellar surface, 
locking away the least abundant element. The remaining 
majority element, either carbon or oxygen, dominates 
the dust chemistry. CO may be dissociated in atmo- 
spher ic shock waves o r by chromospheric UV irradia- 
tion (|Beck et al.l 119921 ). This can lead to the produc- 
tion of small amonuts of other carbon-rich molecules: 
mainly CO2 (obser ved in our stars: see Paper I), but 
also HCN and CS (iDuari et al.lll999D. Small quantities 
of ca rbon-rich dust may then form (|H6fner fc Andersen! 
|2007| ). though oxygen-rich dust production is still ex- 
pected to be dominant. In our sample, however, the 
featureless component often appears to be the dominant 
species, and no evidence for oxygen-rich dust is found. 

Having ruled out the alternatives, we conclude that 
the source of the dust excess must be iron grains, and 
that these must dominate the condensed material in 
some stars. If we naively assume identical populations 
of spherical grains, the fraction of warm silicates in the 
'naked' stars being < 1% of the iron content (Figure |31 
the same would be true at a < 3% level for an amC com- 
ponent). The silicate fraction may increase if iron grains 
are made smaller and oblate, thus presenting a higher 
cross-section per unit mass: extreme needle-like grains 
are required if silicates are to become the dominant dust 
species by mass. As we do not know the grain shape, it 
becomes difficult to estimate a mass-loss rate for an iron 
wind. Spherical grains produce unfeasibly low expansion 
velocities of order ^m s^^ for a pure iron wind; expan- 
sion velocities for NGC 5927 V3 and NGC 6352 V5 were 
not calculated as the model assumes a zero initial veloc- 
ity for both temperature components, which will not be 
the case if they form at different radii. The mass-loss 
rates and expansion velocities can also be significantly 
modified by altering, e.g., the grain size or density (sin- 
gle crystal vs. conglomerate), but reasonable changes 
still fail to produce feasible values, except for extremely 
elongated grains. 

This implies that iron alone may not be able to drive 
a dusty wind in these stars. Other factors may also ac- 
celerate the wind. The se could include small am ounts 
of amorphous carbon (JHofner fc AndersenI l2007t ): en- 
ergy depositio n eithe r from shocks created by stel- 
lar pulsations (lBowen| [T988'). or magnetic reconnection 
(jPiipers fc HearnI 11989*): increased opacity from circum- 
stellar molecular gas (.Elitzur et al.,,1989.) : stellar pulsa- 
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Fig. 3. — Fits to the Spitzer IRS spectrum of uj Cen V42 using alternative emission sources. In each panel, the following is shown: dots — 
IRS spectrum; plus sign — phase-corrected L-band photometry; dotted line — fit for indicated grain composition or emission mechanism. 



tion, either by direct acceleration or dissipation of shocks; 
line-driving of H i or Ca ii atoms, or H2O molecules (e.g. 
lElitzur et all 119891 : lBowenlll989[) . The driving mecha- 
nism is thus likely to be complex, with radiation pres- 
sure on metallic iron being only one of several accelerants 
which drive matter from the star. Obtaining an accurate 
radial profile of the wind's outflow velocity will be crucial 
in determining which factors dominate this acceleration. 

3. DISCUSSION 

3.1. The fate of iron and iron oxide grains 

Efficient production of metallic iron by oxygen- 
rich AGB stars would explain several iron deple- 
tions. Gas-phase de pletion of iron e xtend s to the 
interstellar medium (ISavage fc BohlinI I1979H . Iron 
grains contribute to in frared emission from this medium 
(jChlewicki fc LaureiislligSSi') and iron is s een to be de- 
stroved in shocks (jLebouteiller et aLll2008[ ). Iron is also 
deple ted by > 90% in the gas ph ase of planetary neb- 
ulae (jDelgado Inglada et al.l l2009t ) . Direct evidence for 
iron condensation in AGB winds comes from depletion 
patterns in binary post-AGB stars. In these stars, the 
wind is captured in a circumbinary disk, from where gas 
re-accretes onto the star. These stars show large deple- 
tions of refract ory elements, wit h iron being among the 
most depleted (|Maas et al.ll2005l) thus indicate a history 
of iron-rich dust production. 



The 20-/im feature, seen in NGC 6352 V5, is also 
observed in NGC 5927 VI and Terzan 5 V6. It has 
a width (FWHM) of ^ 3/im and a line-to-continuum 
ratio up to 50%. This feature has previously been 
seen in low mass-loss stars; it correlates with the pres- 
ence of an ll-^m shoulder, and 13.1- and 28-Atin peaks 
(jLittle-Marenin fc Little|[T98l ISloan et al.ll200l . The 
related peaks are also present in our sample, and we 
confirm that the 20-/im feature disappears for more 
evolved (i.e. more luminous) stars. The 20-/im and asso- 
ciated features are barely detectable in the naked stars, 
strongest in stars with weak silicate emission and dis- 
appear when silicates dominate (Figure |4|). This sug- 
gests they are produced either during the formation or 
destruction of silicate grains. The identification of the 
feature with Mga;Fei_a;0 is unconfirmed mainly because 
it does not predict the associated features, but the pro- 
posed identificatio ns for those (corundum, pyroxenes) 
(jSloan et al.l 120031 ) do not provide good fits at 20 ^m, 
making it plausible that a range of minerals form under 
these conditions. We suggest that these features may be 
caused by the chemical modification of silicate grains: sil- 
icate grains which adsorb iron into the lattice have higher 
opacity, and can be heated until the iron-bear ing com- 
ponent dissociates from the grain (jWoitkdl2006f l. This is 
expected to occur at any r a.dius from the star at w hich 
dust formation is efficient (|Gail fc Sedlmavil 119990 . In 
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Fig. 4. — Variation of dust properties with bolometric magnitude 
and nietallicity. Symbols as Figure [l] 

this model, metallic iron is produced first, but at higher 
metallicity and mass loss rates, the iron becomes partly 
incorporated into other grain types including iron-oxide 
and iron-containing silicates. 

Studies of inclusions within meteoric primitive Solar 
System material have traced dust grains to specific ori- 
gins in AGB st ars, supernoyae and novae, through their 
isotopic ratios (|Zinneiil2003l : [Messenger et al.ll2005l ). Re- 
centl y, the first pre-sol ar iron-oxide grain has been identi- 
fied ( Floss et al.ll2008l ). The iron and oxygen isotopic ra- 
tios in it are indicative of an origin in an AGB star. This 
supports the hypothesis that AGB stars are a source of 
interstellar iron dust, and have contributed to the Solar 
System's iron content. 

3.2. The dust condensation sequence at low [Fe/H] 

Figure m shows the variation of dust types (see Paper I) 
with bolometric magnitude and metallicity. Silicate-rich 
winds are only seen at high luminosity or high metallic- 
ity, with iron dominating the outflow in lower-luminosity, 
lower-metallicity stars. This indicates that the dust con- 
dens ation sequence depends on metallicity and luminos- 
ity. IGail fc Sedlmavri (|1999( ) show that, under equilib- 
rium conditions, iron is only expected to condense before 



silicates in regions with high gas pressure. If this applies 
to our stars, it implies high gas pressure at low metallicity 
and/or low luminosity, as may be expected if the outflow 
velocities where iron condenses (at the outside edge of the 
molsphere) are very low. Under non-equilibrium condi- 
tions, gr ain growth chem i stry i s strongly controlled by 
kinetics (jTachibana et al.ll20Q9[ ). which may favour the 
for mation of iron. 

In lGail fc Sedlmavil ()1999f) , silicates form before iron at 
lower densities and temperatures. The restriction of sil- 
icates towards higher metallicity and luminosity (Figure 
[3]) suggests that their circumstellar envelopes are more 
extended, favouring silicates. Silicate production is ex- 
pected to increase with b oth stellar metallicity and lu- 
minosity (|van Loonll200ClD . 



4. CONCLUSIONS 

Using mid-infrared spectra and photometry, we 
have shown that there is a considerable amount of 
unattributed infrared flux in a large selection of globu- 
lar cluster giants. Through radiative transfer modelling, 
we deduce that this flux is most likely due to metallic 
iron grains forming in a truncated stellar wind, corrobo- 
rated by a potential identification of FeO with the 20-/im 
emission feature. The production of metallic iron seems 
to become more dominant at lower metallicity, suggest- 
ing that the dust condensation sequence is fundamentally 
different in metal-poor objects. Large-scale production 
of iron grains and iron oxide in AGB stars can explain 
iron depletion in the gas and solid phases of the post- 
AGB stars and planetary nebulae; as well as isotopic 
ratios in FeO grains in meteorites. While iron increases 
opacity in oxygen-rich winds, it remains unclear whether 
it can dominate driving of metal-poor winds. Future ob- 
servations to determine outfiow velocities in these stars 
should help determine the driver of these stellar winds. 

This paper uses observations made using the Spitzer 
Space Telescope, operated by JPL, California Institute 
of Technology under NASA contract 1407 and supported 
by NASA through JPL (contract number 1257184). We 
thank Martha Boyer for her helpful comments. 
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